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The large J2 limit of the square-lattice Ji — J2 Heisenberg antiferromagnet is a classic example of 
order by disorder where quantum fluctuations select a collinear ground state. Here, we use series 
expansion methods and a meanfield spin-wave theory to study the excitation spectra in this phase 
and look for a finite temperature Ising-like transition, corresponding to a broken symmetry of the 
square-lattice, as first proposed by Chandra et al. (Phys. Rev. Lett. 64, 88 (f990)). We find that 
the spectra reveal the symmetries of the ordered phase. However, we do not find any evidence for a 
finite temperature phase transition. Based on an effective field theory we argue that the Ising-like 
transition occurs only at zero temperature. 
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The square-lattice J\ — Ji Heisenberg antiferromagnet 
is described by the Hamiltonian 

n.n. n . n . n . 

where the first sum runs over the nearest neighbor and 
the second over the second nearest neighbor spin pairs of 
the square-lattice. For J\ — 0, the two sublattices are dis- 
connected and individually order antiferromagnetically. 
For non-zero J1/J2 (less than some critical value), in the 
classical ground state, the two sublattices remain free to 
rotate with respect to each other. However, quantum 
fluctuations lift this degeneracy and select a collinear or- 
dered state, where the neighboring spins align ferromag- 
netically along one axis of the square-lattice and antifer- 
romagnetically along the other [1-4]. Thus the ground 
state breaks both spin-rotational symmetry as well as 
the four-fold symmetry of the square-lattice. 

It is well known from the Mermin- Wagner theorem 
that the continuous spin-rotational symmetry is restored 
in 2D at any finite temperature. However, the discrete 
broken lattice symmetry can in principle survive at finite 
temperatures. In an important paper, Chandra, Cole- 
man and Larkin [2] argued that this symmetry should be 
restored at a finite temperature phase transition, which 
lies in the 2D Ising universality class, and gave estimates 
for the transition temperature. They made use of a lin- 
ear spin- wave theory (LSWT) for the excitation spectra, 
which had gapless excitations at four points of the Bril- 
louin zone ((0,0), (n,0), (0,7r) and 

Recent interest in this model comes from the discovery 
of two materials L^VOSiC^and L^VOGeO^y Melzi et 
al. [5,6]. Electronic structure calculations for these mate- 
rials [7,8] lead to J2 bigger than Ji, perhaps by as much 
as an order of magnitude. The reason for the unusually 
large J2 can be qualitatively understood from the crystal 
structure. The VO5 pyramids alternately point up and 
down, and hence, the spin-half vanadium atoms are al- 
ternately displaced slightly above and slightly below the 



plane. This causes an increased overlap with the second 
neighbors which fall in the same plane. Many experi- 
mental features of the material are well described by the 
Ji — J2 Heisenberg model. However, the value of the ratio 
J1/J2 remains ill-determined [9]. 

It was argued by Rosner et al. [8] that the measure- 
ment of the spin-wave spectra maybe particularly useful 
for establishing this ratio. One of the purposes of this pa- 
per is to present quantitatively accurate spectra for the 
model going beyond linear spin-wave theory. Indeed, as 
we will show, the full spectra have gapless excitations at 
only two symmetry related points of the Brillouin zone 
(0,0) and (0,tt), whereas the accidental degeneracies at 
(tt, 7t) and (ir, 0) are lifted by the order by disorder effect 
[1] (here we choose the x axis along the direction of the 
collinear spin ordering, see Fig. 1). Furthermore, the 
ratio of the spin- wave velocities along x and y directions 
depends sharply on the J1/J2 ratio and its measurement 
can indeed be used to determine the latter ratio. 

The real materials also have weak interplanar cou- 
plings, which eventually lead to 3D long-range order and 
a finite temperature transition. The specific heat data 
above the 3D transition does not have any clear features 
which could be interpreted as a 2D Ising transition. In- 
deed, the issue of finite-temperature Ising-like transitions 
may only be relevant to materials with sufficiently weak 
interplanar couplings. The issue is conceptually impor- 
tant, however. Can long-range Ising order, corresponding 
to antiferromagnetic correlations along one axis and fer- 
romagnetic correlations along the other, survive at finite 
temperatures even after the underlying spin-spin correla- 
tions become short ranged? That is the somewhat para- 
doxical prediction of Chandra et al. [2] 

We calculate the susceptibility appropriate for this 
Ising order by high temperature series expansions. The 
numerical series extrapolations fail to substantiate a fi- 
nite temperature Ising-like transition, they are only con- 
sistent with a T = phase-transition in a strictly 2D 
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model. We also present field-theoretical arguments in 
favor of a T — phase-transition, and find the critical 
index for the T = divergence of the susceptibility. 

We begin by calculations of the excitation spectra at 
zero temperature. The standard linear spin-wave theory 
predicts the dispersion relation [2], 
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and this results in a spectrum with zero modes at four 
points k = (0, 0), (0, tt), (tt, tt) and (tt, 0), as discussed by 
Chandra et al. [2]. 

We have derived the spin-wave dispersion using a 
meanfield spin- wave theory (MFSWT), see e.g. Refs. 
[4,10]. The result reads 
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where fj,, x and y are functions of J^/ J\ . For example for 

J 2 /Ji = 2 : n = 1.153, x = Q.327, y = Q.174 
J 2 /Ji = l: M = 1.122, a: = 0.783, y = 0.452 (4) 
J2/J1 = 0.6 : fi = 0.818, x = 1.583, y = 1.323 

This dispersion give a spectrum with zero modes at only 
two points, i.e. k = (0,0) and (0,tt). 

We have performed an Ising expansion [11] of the ex- 
citation spectra using the linked-cluster series expansion 
method, as reviewed recently by Gelfand and Singh [12]. 
The spin triplet excitation energy A(k x ,k y ) has been 
computed up to order 9 for the system in the collinear or- 
dered phase. A list of 1796200 linked clusters of up to 10 
sites contribute to the triplet excitation spectrum. The 
series are available on request. Note that the dispersion 
relation has the symmetry: A(k x , k y ) = A(k x , tt — k y ). 

The excitation spectra for coupling ratios J2/J1 = 
2,1,0.6 are shown in Fig. 1, together with those ob- 
tained from standard linear spin-wave theory [2] and a 
MFSWT [10]. We can see from this figure that the spec- 
trum does not have the 90° rotational symmetry of the 
square-lattice. There are zero modes at only two points 
i.e. k = (0,0) and (0,tt), as predicted by the MFSWT 
theory, but not the standard theory. In fact, the series 
expansion data agree with the predictions of the MFSWT 
theory with remarkable accuracy over the entire range of 
momenta, except for a small dip in the third panel. The 
gaps at k = (tt, tt) and (tt, 0) are non-zero, due to the 
"order by disorder" effect [1]. 

The gap at k = (tt, tt) and the spin-wave velocity along 
x and y-directions versus J1/J2 are shown in Fig. 2 and 
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FIG. 1. Plot of the spin-triplet excitation spectrum 
A(k x ,k y )/ Ji along high-symmetry cuts through the Brillouin 
zone for the system with coupling ratios J2/J1 = 2,1,0.6 
shown in the figure from the top to the bottom, respectively; 
the blue dotted lines are the results of MFSWT, the green 
dashed lines are the results of LSWT, while the red points 
with error bar are the results of series expansions. The un- 
perturbed spin configuration assumed in x and y directions 
are shown as inset. 



Fig. 3. Clearly, the latter anisotropy provides a good 
way to determine the J1/J2 ratio. 

We now turn to the question of an Ising- like finite tem- 
perature transition. To explore the possibility of such a 
transition we calculate the high temperature series for 
the susceptibility with respect to the field 
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i.e. we compute the series for 

T X = (F 2 ) - (F) 2 



where 
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and P = l/(k B T); we take k B = 1 here. Note that (F) 
is zero for the bulk system, but we need to include it in 
Eq. (6) because it is not zero for each individual cluster 
in the series expansion. 

This field breaks the 90 degree rotational symmetry, as 
in the zero temperature collinear ordered phase, and if 
there was a finite temperature transition, one would ex- 
pect x t° diverge at the critical temperature. The series 
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FIG. 2. The spin-triplet excitation gap A/J2 at k = (n, n) 
versus Ji / J2 obtained from series expansion (the points) , and 
the MFSWT (the solid line). 
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FIG. 3. The spin- wave velocity V/J2 along x and 
y-directions versus J1/J2 obtained from series expansion (the 
points with dotted lines connecting them), and the MFSWT 
(the dashed line). 



has been computed up to order /3 8 , using the same list of 
1796200 linked clusters of up to 10 sites as for the spin 
dispersion. The series are available on request. 

For the full series, we first tried to locate the critical 
point by Dlog Pade approximants [13], but this did not 
give any consistent results. We also used integrated dif- 
ferential approximants [13] to extrapolate the series, and 
the results are shown in Fig. 4. We can see that 1/x 
appears to vanish only at zero temperature, rather than 
at finite temperature. This is consistent also with the 
results of the Pade approximants. The Ising-like transi- 
tion temperatures estimated according to [2] are shown 
in Fig. 4 by arrows. 

To look further into this issue, we have also re-analyzed 
our high-temperature series for the specific heat C v [7,8] , 
for the ratio of couplings in the collinear ordered phase. 
While the convergence of the series extrapolation is not 
very good at low temperatures, the Dlog Pade approx- 
imants to the series for C v do not show evidence for a 
finite temperature critical point [9]. 

We also analyse our numerical results using an effec- 
tive field theory. The dynamical variables of the problem 
are two vector fields ni and n 2 corresponding to two 
interpenetrating Neel sublattices with exchange interac- 
tion J 2 . Each of the fields is described by the nonlinear 
er-model, and there is also an interaction between the 
fields. The effective Lagrangian of the system reads [2] 
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where g ~ J\ /J 2 is the coupling constant. There arc two 
degenerate energy minima corresponding to a = ni • n 2 = 



±1. So a imitates an Ising parameter, and at T = it 
takes a values corresponding to one of the minima, a = 1 
(x-collinear state) or a = — 1 (y-collinear state). How- 
ever, we stress that there is not a dynamical Ising vari- 
able in the system, the entire dynamics being described 
in terms of ri! and n 2 and the Lagrangian (8). The 
two minima a = ±1 are separated by a potential bar- 
rier W(Q) <~ EL 2 , where E ~ g is the interaction energy 
density, and L — > 00 is the size of the system. 

At a finite temperature, each sublattice behaves as 
an ordered Neel state up to a length scale £ sa 
0.303 exp(27rp s /T)[l - T/(4wp s )}, where p s « 0.182J 2 is 
the spin stiffness [14]. At low temperature the energy 
barrier required to move from the a = +1 configura- 
tion to the a = — 1 configuration, W(T) ~ E£ 2 , remains 
finite at any finite T! So, using the quantum field the- 
ory language, one has to say that there is a nontrivial 
instanton in the problem with tunneling probability oc 
cxp(— E£ 2 /T). This is different from the true Ising situ- 
ation, where there is a dynamical Ising field with stiffness, 
and the tunneling probability oc cx.p(-E(T)L 2 /T) — > 0, 
as L — > 00. The "true" Ising transition occurs when E{T) 
vanishes. 

Based on the picture described we immediately come to 
the conclusion that at any finite T the system fluctuates 
beween a — 1 and a = — 1 local minima, and hence 
the transition temperature for the Ising-like transition 
is zero in agreement with numerical data. According to 
this picture one would expect that at low temperature 
and at J 2 > J\ the suceptibilty with respect to the field 
(5) behaves as 
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FIG. 4. The inverse of susceptibility 1/JiX versus T/Ji for 
J1/J2 ~ 1,0.25,0.1. Several different integrated differential 
approximants to the high-temperature series are shown. The 
fsing-like transition temperatures estimated according to [2] 
are shown by arrows. 



where a is a critical index, (we neglect the prefactor 
[1 - T/(4irp s )} in £). Plots of ln(Jix) versus J 2 /T are 
presented in Fig. 5. The plots are roughly consistent 
with prediction (9) with the critical index awl. 

In conclusion, in this paper we have studied the 
collincar phase of the J\ — J2 square-lattice Heisen- 
berg antiferromagnet by series expansion methods. We 
have obtained quantitatively accurate excitation spectra, 
which are significantly more accurate than the spin- wave 
calculations [2], and should be helpful in determining 
the exchange parameters for materials with large second 
neighbor interactions. The MFSWT predicts these spec- 
tra with great accuracy. We have also studied an Ising- 
like phase transition in this model. Using high tempera- 
ture expansions for the appropriate susceptibility as well 
as quntum field theory arguments we have shown that 
the transition temperature is exactly zero. 
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FIG. 5. In(Jix) versus J 2 /T for J x /J 2 = 1,0.25,0.1. Also 
shown are 0.8 + ln£ with (red bold line) and without (blue 
straight line) the prefactor [1 — T /{Air p s )]. 
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